From the self-force problem to the Radiation reaction formula 



in 



Yasushi Mino * 

mail code 130-33 Caltech Pasadena CA 91125 USA 
February 7, 2008 



Abstract 



We review a recent theoretical progress in the so-called self-force problem of a general relativistic two-body system. 
Although a two-body system in Newtonian gravity is a very simple problem, some fundamental issues are involved 
in relativistic gravity. Besides, because of recent projects for gravitational wave detection, it comes to be possible to 
see those phenomena directly via gravitational waves, and the self- force problem becomes one of urgent and highly- 

■ motivated problems in general relativity. Roughly speaking, there are two approaches to investigate this problem; the 
^ so-called post-Newtonian approximation, and a black hole perturbation. 

C\| ■ In this paper, we review a theoretical progress in the self-force problem using a black hole perturbation. Although 
the self-force problem seems to be just a problem to calculate a self-force, we discuss that the real problem is to define 

■ a gauge invariant concept of a motion in a gauge dependent metric perturbation. 

T— I ■ I. INTRODUCTION 

m ; 

I ' For a single isolated body, the self- induced gravitational field does not accelerate the body because of the momentum 
^ . conservation. However, in a multi-body system in general relativity, there could be the case that the self-induced 
(~\| ■ gravitational field could accelerate the body. Although one may naively say we have the self-force acting on the body, 
, the concept of the body's motion is itself ambiguous in general relativity because one can freely choose a reference 

■ coordinate system to describe it. Thus, even if we could derive an explicit form of the self-force, there still remains 
' a fundamental question; what we learn from the the self-force about the body's motion. Here we focus on a binary 

■ system as a simple example to consider this fundamental question. 

[ The dynamics of a binary system is also motivated by rapidly-progressing projects of gravitational wave detection. 
— The so-called post-Newtonian approximation has a fairly long history in calculating an evolution of a equal-mass 
^ I binary system. Although we now have a huge amount of interesting results, the post-Newtonian approximation was 
I usually done only with a few choices of the coordinate system. 

Because the choice of coordinate systems is a key issue in understanding this fundamental question, it would be 
more interesting to consider another approach in which we may have a more degree of freedom to choose coordinate 
K>" ■ systems. Here we consider the self-force problem by a black hole perturbation. We consider the relativistic two 
vTji I body problem as a background black hole and a point particle orbiting around the black hole. We assume that the 
, background black hole is described by a Kerr geometry, and consider a metric perturbation induced by a particle 
' moving along a general bound orbit. The reference coordinate system for the particle's motion is dominantly given 
by the coordinate system of the background Kerr geometry, but, since the background geometry is perturbed by the 
self-gravity of the particle, we have a small freedom to choose the coordinate system. This is what we know as a 
gauge freedom of the metric perturbation. 

This two-body system is considered to be a good approximation of a gravitational-wave source we may see by the 
LISA project. Because of the long base line, LISA will see low frequency gravitational waves (around lOmHz to O.lHz). 
Gravitational waves in this frequency range may be obtained from an orbit around a supermassive black hole (around 
10^ to 10* Mq), which we may astrophysically expect to exist at the center of a galaxy. Our recent understanding 
of a dynamical structure of a galactic core suggests that supermassive black holes in the observable range of LISA 
would capture stellar mass compact objects in their gravitational potentials. Those gravitationally captured objects 
are inspiralling around the black holes, emitting gravitational waves detectable by LISA. 
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II. OVERVIEW 



The study of an orbital motion of a particle in a black hole geometry has started from the discovery of the so-called 
MiSaTaQuWa self- force [1], which wc shall review in Sec. III. At this moment, the question was how we deal with a 
point particle and we obtained a general covariant expression of the self-force in a given vacuum background. Since 
the MiSaTaQuWa self-force is just a formal result, our next question was to evaluate the self-force explicitly in a Kerr 
geometry, using a known technique to calculate a metric perturbation. We review some basic ideas of some promising 
calculation methods in Sec. IV. 

There are two approaches. One is to calculate the self-force directly by explicitly implementing the point splitting 
regularization of the MiSaTaQuWa self-force, and we call it the self-force calculation. The other is to calculate 
the radiation reaction to the constants of motion based on our understanding of the self-force, and we call it the 
radiation reaction formula. By the self-force calculation, one may have the exact self-force in the framework of a 
linear perturbation, however, we still have a technical problem when we consider the self-force calculation in a Kerr 
black hole. Besides, as we shall discuss in Sec.V, there is a theoretical issue whether the self-force derived by this 
method correctly describes the radiation reaction effect. On the other hand, the radiation reaction formula considers 
radiation reaction to the orbit. Different from the self-force calculation, we only calculate the radiative part of the 
self- force in a class of gauge conditions where we are guaranteed to include the radiation reaction effect to the orbit. 
It has a technical advantage that the calculation method of this formula is already established. Furthermore, our 
recent result shows that the radiation reaction formula may have the most optimal prediction of the orbit. 

It is widely believed that the self-force describes the orbital evolution with the effect of gravitational radiation 
reaction. Sec.V discusses that this is not true in general. We point out that a usual scheme of a metric perturbation 
describes a binary evolution only for a finite time interval. We show that the self-force could vanish by a special 
gauge transformation for this whole time interval, irrespective to gravitational radiation we see from the asymptotic 
gravitational wave flux. Thus, strictly speaking, the self-force has nothing to do with gravitational radiation reaction 
in the usual metric perturbation without some restriction of gauge conditions. We consider this problem happens 
not because of our definition of the self-force but because the usual scheme of a metric perturbation cannot describe 
the binary evolution long enough. In order to solve this, we propose an adiabatic extension of a metric perturbation, 
which we review in Sec. VI. The adiabatic extension is possible under a class of gauge conditions, where the self-force 
is guaranteed to include the radiation reaction effect. We also discuss the validity of the adiabatic approximation of 
the linear metric perturbation and the self-force defined by this. We find that the self-force in this class of gauge 
conditions might not predict the binary evolution sufficiently long for the LISA project in general. In order to extend 
the validity, we introduce a special gauge condition where the radiation reaction formula gives an exact self-force. In 
this gauge condition, we find that the self- force may predict the binary evolution long enough for the LISA project. 

III. SELF-FORCE 

In this section, we review a derivation of the MiSaTaQuWa self-force. In this review, we focus on the physical 
motivation why we take our calculation strategy and we omit technical details of the calculation. We suggest the 
readers interested in those details to refer our original paper, or, the review papers by us or Poisson [1]. 

A. A historical motivation 

In a black hole perturbation, a point particle is not well-defined because the self-gravity of the point particle 

diverges and the metric perturbation becomes invalid around the particle. Nevertheless, wc had been using the point 
particle without any rigorous supporting argument because of the technical simplicity and because the point particle 
might be an only possible description of a star whose internal structure is supposed to be unimportant, and, a metric 
perturbation induced by a point particle has been investigated for a long time [2] before the MiSaTaQuWa self-force 
was derived [1]. 

At this time, we used the so-called Press- Teukolsky formalism [3] which gives us a finite energy and angular 
momentum flux carried away by gravitational waves. Even if the point particle induces the divergent gravitational 
perturbation, one can define a finite energy and angular momentum on a t-constant hypersurface in the background 
Kerr geometry and we consider that the loss of the energy and angular momentum by the gravitational wave emission 
describes radiation reaction to the particle's motion. This method is now refered by the balance formula. 

The balance formula was considered to be successful to describe radiation reaction to a particle moving along a 
quasi-circular orbit or an equatorial orbit. In this case, we consider an adiabatic approximation of the orbit where 
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the orbit is approximated by a geodesic at each instant. Then, instead of the evolution of the orbital coordinates, we 
consider the evolution of "3 constants of motion" (the orbital energy E, z component of the angular momentum L and 
Carter constant C) of geodesies around a Kerr geometry. The balance formula tells us radiation reaction to only two 
of these constants, E and L, and the problem was to derive radiation reaction to the rest, C. When the orbit is either 
quasi-circular or equatorial, one can prove that it will remain quasi-circular or equatorial under radiation reaction [4] , 
as a result, we have an algebraic relation with E, L and C to derive radiation reaction to C by the balance formula. 

The difficulty to compute radiation reaction to C was that Carter constant is not associated with the Killing 
vector. The Press- Teukolsky formalism was derived because E and L are the conserved currents defined on a t- 
constant hypersurfaces in association with the symmetry of the background spacctime. On the other hand, Carter 
constant is defined only along a geodesic and has nothing to do with the symmetry. Besides, it came to be aware that, 
in order to integrate the geodesic equation, one needs three phase constants in addition to the constants of motion, 
and we did not know how those phase constants evolve by gravitational radiation reaction. 

At this moment, a simplest solution of this problem seemed to be an explicit calculation of the self-force. We 
considered that the self-force must have contained all the information of the orbit with the effect of radiation reaction, 
and that, after computing the orbit using the self-force, one can calculate a gravitational waveform by a linear 
metric perturbation formalism. As we shall argue in Sec.V, this statement is proven to be incorrect because a usual 
perturbation formalism cannot describe radiation reaction to the orbit, nevertheless, this idea had been spread widely 
in the gravitational wave community. 



B. Problem of Divergence 

A point particle is defined by the Dirac's delta fimction. It has no volume, but, a finite mass. This makes a 
divergent metric perturbation at the location of the particle. Because of the divergence, we have two problems; 

1) Can we use the point particle to induce a linear metric perturbation, which invalidates the linear perturbation 

scheme? 

2) Can we extract a finite physical information of a self-force from the divergent field? 

A common idea to understand a divergent field in Quantum Field Theory is renormalization. In electromagnetism, 
a point charge induces a divergent vector potential at the location of the particle. When the charge moves in a curved 

spacctime, this self-induced vector potential will act on the charge's motion. Because the quantum electrodynamics is 
renormalizable, a formal derivation of this electromagnetic self- force in a curved spacctime was successfully made by 
a mass renormalization [5] . Here one evaluates a surface integration of the stress-energy tensor of the charge and the 
electromagnetic field over a world tube surrounding the orbit. The sum of the surface integral must vanish because 
of the gauss's law, and we obtain the self-force equation including the effect of electromagnetic radiation reaction. 

We considered an analogous calculation for a gravitational self-force by the mass renormalization [1]. We do not 
have the stress-c^ncrgy tensor of the particle and the full gravitational field, however, we could find a conserved current 
of the particle and the metric perturbation in the background geometry. We consider to expand the Einstein tensor 
by the metric perturbation h fj,!/ in the background. Qui/ as 

G^. [g + h] = G,^ [g] + [h] + [h] , (3.1) 

where gI/J and G^^i^^ represent the linear terms and all the non-linear terms of the expansion with respect to /i^^. 

Because we assume the background is a vacuum solution, we have Gp,i,[g] = 0. It is well-known that the divergence 

of G^,n, [h] with respect the background geometry vanishes algebraically. Thus, we have the conserved current in the 

background geometry as 

V = G{f+)[/i]+T^,, (3.2) 

where T^i, is the stress-energy tensor of the particle. Because G^^i^^ [h] is quadratic in /i^^ to the leading order of 
the metric perturbation, we may consider this as an "effective stress-energy tensor" of the particle and the metric 
perturbation defined in the background geometry. By doing a surface integration over a world tube, we obtained a 
gravitational self-force equation by the mass renormalization. 

Although we could have derived the MiSaTaQuWa self-force by the mass renormalization, the derivation did not 
answer the above fundamental questions. A crucial issue in this derivation is that the perturbative quantum gravity 
is known nonrenormalizable, thus, there is no theoretical support to use the mass renormalization in deriving the 
gravitational self-force. Furthermore, a physical interpretation of the motion is not clear yet. If we suppose that the 
"point particle" represents a black hole, for example, the center of the particle is inside the event horizon, and one 
may ask what it means to consider a motion hidden inside the horizon. 
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C. Matched Asymptotic Expansion 



In order to answer these fundamental questions on the use of the "point particle", we consider another approach. 
Since the primary problem here is a theoretical justification to use the "point particle" rather than just a calculation 
of the self-force, we consider to construct a metric itself by using a matched asymptotic expansion technique. By 
the matched asymptotic expansion, we consider to locate a small black hole of a length scale m in a background 
spacctimc of a curvature scale L. Assuming that m << L, one could construct the metric perturbativcly by a small 
parameter m/L << 1. This approach was originally taken by the post-Newtonian calculation [6], where we have an 
explicit expression of a metric easy to control. In this problem, we consider to locate the small black hole in a general 
vacuum background spacetime and, in order to have an explicit expression of the metric, we use the so-called bi-tensor 
formalism introduced in Ref. [5]. 

With a radial distance r(> 0) from the center of the small black hole, the gravitational field is determined dominantly 
by that of the small black hole at r < L, and by the background spacetime at r > m. Hence, the metric at r < L is 
constructed by the small black hole with a black hole perturbation (we call the near-zone expansion) as 

= + L-'al^J' + L--g^^P +■■■, (3.3) 

where g^^ is the metric of the small black hole. The black hole perturbation (/^"^^ is defined with an appropriate 

boundary condition on the black hole horizon. As we cannot specify a boundary condition for glZ/"^ at r > L, 
the metric perturbation cannot be determined uniquely until we complete the matching. The metric at r > m is 
constructed by the background spacetime with its perturbation (we call the far-zone expansion) as 

= 9^1 + wffi/Jl + m2<7l/Jl + • • • , (3.4) 

where jg the metric of the background spacetime. The perturbation of the background gllu^ is derived with 
an appropriate asymptotic boundary condition and is not determined uniquely because of the lack of a boundary 
condition at r < m. Instead of the boundary condition at r < m, one can put an auxiliary matter at r < m to induce 
a vacuum perturbation at r > w, and we use a point source with an arbitrary internal structure derived in Ref. [7] . 

The matched asymptotic expansion can be done at the radius where both near-zone and far-zone e xpan sions are 
valid, which we call the overlapping zone. Suppose we take the radius of the overlapping zone as r VmL, and one 
can re-expand (3.3) and (3.4) as 

<r = (9^^ + mg%^'^ + ...) + L-\glT^ + m^,!'^ + ■■■) + ■■■, (3.5) 

ff^r = + +•••)+ m{g^T^ + L-^gll'^ +■■■) + ■■■. (3.6) 

(3.7) 

From a dimensional analysis, we may see g^^^J"^^ ^ ffjli/'"' ~ r""™, thus, at the overlapping zone, one can estimate 
the expansions as 

I,-"m™5("/'") ~ L-^m""^!"/""! ~ (m/L)("+'")/2 . (3.8) 

Under the matched asymptotic expansion, we assume mjL << 1 and we equate g^^^J™^ = g^i^J"^^ one by one from 
small to large n + m which will perturbativcly determine the boundary condition of (3.3) at r > L and the internal 
structure of the auxiliary matter to induce the perturbation of (3.4). 

Here, it is important to note coordinate conditions. Since the derivation of (3.3) and (3.4) are totally independent, 
they may use different coordinate systems, and we have to introduce a coordinate transformation so that we can 
correctly evaluate g^^J™^ = 5}!]/'"' . It is also important to fix a gauge condition for the metric of the small black hole 
in order to fix the location of the particle in the near-zone expansion. 

We suppose that the background metric and its coordinate system are given by {g^^, {2^"}) and that the orbit in 
the background metric is described by 2:"(t) with a parameter r. We also suppose that the black hole metric and 
its coordinate system are given by {gfj^,{X^,i = 0, 1,2,3}) and we take a Schwarzschild black hole metric in the 
harmonic coordinates where X° is the temporal coordinate and X^,X'^,X^ are spatial. In order to fix the location 
of the small black hole in the near-zone expansion, we take the gauge condition such that the I = 1 modes of the 
perturbation L~"g^^^^ in the tensor-spherical decomposition vanish. Then, we may say the particle is located at 
X^ = X^ = X^ = in the near-zone expansion. 
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For the matched asymptotic expansion of the metrices in the far-zone expansion and the near-zone expansion, we 
define the coordinate transformation by using a tetrad defined along the orbit and derivatives of the bi-scalar of 
the half-the-squared geodesic distance. The coordinate transformation is determined perturbatively by the matched 
asymptotic expansion of the metrices. 



D. MiSaTaQuWa Self-force 

We find that the coordinate transformation gives us the information of the particle's orbit because it relates the 
location of the small black hole in the near-zone expansion (i.e. — X"^ = X^ = 0) to the coordinates in the far-zone 

expansion. By the matching g'fj^'' = ^j?/"^ and ^ = g)ll^^ , we find that the particle moves along a geodesic. By 

the matching g'^^J^^ = gfl°\ g^^J^^ = g^^l^\ glL^^^ = ff^V^' and gj^J^^ = gfl^\ we derive the self-force. 

The point particle is shown to be consistent for the derivation of the metric in the far-zone expansion as a result of the 
matched asymptotic expansion of the metric. Even though the point particle makes a divergent metric perturbation, 
it is just an auxiliary source to induce a finite metric perturbation at r > m in the far-zone expansion. This justifies 
the use of a point particle in calculating an energy and angular momentum flux via an asymptotic metric perturbation 
since it uses the metric perturbation only in the far-zone expansion. 

We also consider that the derivation of the self-force by the matched asymptotic expansion suggests a physical 
interpretation of the particle's motion. In the end, the motion of the particle is considered by the coordinate system 
in the far-zone expansion which is the background metric plus a small perturbation induced by the small black hole. 
For an observer far from the particle, the background metric dominantly determines a reference frame to measure the 
motion of the particle with, whatever non-perturbative things happen around the small black hole, thus, the self-force 
derived by the matched asymptotic expansion is expected to be physically meaningful for such an observer. Because 
the small black hole is small compared to the curvature length of the background metric, one can assign the location 
of the particle in the coordinate system of the regular background metric, and consider the motion of the particle 
without divergence. Here we have a small ambiguity in this assignment of the particle's location because of a finite 
size of the small black hole, which can be interpreted as a gauge ambiguity of the self-force. The location of the 
black hole in the near-zone metric is determined by the gauge condition of the metric perturbation L~^''g'^JK In the 
present derivation of the self-force, the gauge conditions of L~^g^^/^^ and L~'^g^^^'' are important. By the matched 
asymptotic expansion, they fix the gauge condition of the linear metric perturbation mgHg of the far-zone expansion. 



IV. REGULARIZATION CALCULATION 

Since the MiSaTaQuWa self-force is just a formal expression of a self-force in an arbitrary background, the next 
problem was to explicitly calculate a self-force in a specific background, especially, a Kerr black hole for an application 
to LISA project. The formal result of the self-force [1] can be written as 

F"(t)= lim F"[h''^B(^x)], (4.1) 

x—>-z{r) 

F«[/c(a;)] = dT'§^[k{x)]v'^v'^ + v''6Ti3js[k{x)yv^v^ , (4.2) 

C/(^) = MW-/i:7(^)' i^"^^"ir)), (4.3) 

where is a field point, z'^(t) is an orbit, and v'^ is the 4- velocity. ST^^ is a linear tensor differential operator acting 
on a metric perturbation to calculate a perturbed Christoffel symbol. ha0{x) is a linear metric perturbation induced 
by a point particle, and /i^™^(a;) is a singular part derived by the Hadamard expansion. /i^^(a;) is called a regularized 
metric perturbation. Both hap{x) and /i^™^(a;) arc divergent along the orbit, but, ^^'^(•t) is regular along the orbit 
since the divergence of hafjix) and h^a0^{x) cancel each others. Thus, the calculation of the self- force can be done by 
first evaluating it at a field point x", then taking the limit x" — » z"{t). By doing this point splitting regularization, 
we have a finite result. 

There are two promising approaches. One is to calculate the self-force directly by doing a Harmonic mode de- 
composition, and we call it the self-force calculation. The other is to calculate radiation reaction to the constants of 
motion based on our theoretical understanding of the self-force, and we call it the radiation reaction formula. 
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A. Self-Force Calculation 



Here we briefly review a method to compute a self- force directly. The basic idea here is that, since it is difBcult to 
compute F°'[h^^s{x)] of (4.1) directly, we consider to calculate F"[h{x)] and F" (a;)] , then, the self-force can be 
derived by 

F"{t)= \im {F"[h{x)]-F"[h'''"<'{x)]). (4.4) 

x—*-z{r) 

We call F°'[h{x)] and F"[/i^™^(x)] by a bare self-force field and a singular self-force field, which is defined at a field 
point not on the orbit z"(t). It is important to note that the subtraction of (4.4) must be done before taking the 
coincidence limit z°'{t) since F°'[h{x)] and F°'[h^^"^ (x)] are divergent along the orbit. 

This type of calculation is called regularization. Various techniques for regularization calculations are well-known 

in quantum field theory, but we cannot use those techniques bccaiisc of two reasons; 

1) In a usual regularization of quantum field theory, we consider a theory of a global Lorentz invariance and we can 
expand the divergent field by the Fourier mode, respecting the Lorentz symmetry. However, a black hole spacetime 
does not have such a symmetry and we may only decompose hafj{x) into Fourier- Harmonic series. 

2) h"^^^^ {x) is derived by the Hadamard expansion, thus, we have its local divergent expression only around the orbit. 
Besides, in the case of a Kerr black hole as a background, a method to calculate a linear metric perturbation was 
proposed just recently in Ref. [8], and it has not been coded yet. Another simple technique to calculate a metric 
perturbation was proposed in Ref. [9], however, it is known to be valid only for a homogeneous metric perturbation and 
it cannot describe a metric perturbation around the orbit. It is also important to note that the original self-force was 
derived in the harmonic gauge condition, thus, we may have a residual divergence by a difi^erence of gauge conditions 
unless we calculate haf}{x) in the harmonic gauge condition. 

The most successful method at present is the so-called mode; decomposition regularization [10], which is actually 
based on two ideas. One is to use a local coordinate expansion originally proposed by us [11], and the other is to use 
the so-called regularization parameters originally proposed by Barack, Burko and Ori [12]. Here we briefly review a 
basic idea of this regularization calculation in the case that the background is a Schwarzschild black hole. There is 
also some progress on an extension to a Kerr black hole and a method to deal with the residual divergence by a gauge 
condition, but, we do not discuss it here because of its technicality. 

We use the Schwarzschild coordinates and denote a field point and an orbital point by x" = {t, r, 6, (j)} and 
z" — {to, '''O: ^0; '/'o}- We suppose that a full metric perturbation is derived by a harmonic decomposition and we 
calculate a harmonic decomposition of the bare self-force field. ^ 

After the harmonic decomposition, the expansion coefiicients become finite functions of (t, r), and the divergence 
of the bare self- force field appears only by taking the infinite sum over harmonics at the coincidence limit a;" — > 2:". 
We take the angidar coincidence limit [d, (f) {6q, (po), and take the sum over the azimuthal modes. Then the mode 
decomposition of the bare self-force field can schematically be written as 

F"[h]=^F[''^''^^"{t,r). (4.5) 

We note that the harmonic decomposition is ill-defined on a 2-sphere of {t,r) = {to,ro) because the bare self- 
force field diverges, as a result, the mode decomposition is not continuously defined across this 2-sphere as 

As for the singular self-force field to be subtracted, we only have its local expression by the Hadamard expansion, 
and it is not defined on the entire 2-sphere {6, (j>), hence, an exact harmonic decomposition cannot be defined. However, 
the regularization calculation (4.1) depends only on the local behavior of the fields around the orbit, and one can 
freely add a finite term to the singular self- force field if it vanishes at the coincidence limit x" ^ z". We consider the 
local coordinate expansion of the singular self-force field. We suppose that x" — z" is sufficiently small and expand 
the singular self- force field by the power of x" — z" as 



^ Since we are interested in calculating a self- force vector, there are some options in the harmonic decomposition. In order to 
avoid a complication of multiple components, we consider it easier to decompose each component of the self-force vector by the 
scalar harmonics. For this technical detail, see Ref. [10]. 
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= ^ /(« ,„„)^(t - torir - rone - 9or{ct> - <t>or* , (4.6) 

e = IgafsizXx'' - z'^){x^ - z^) , (4.7) 

where gaffiz) is the background metric evaluated at the orbital location. For the regularization calculation, it is 

sufRcicnt to decompose the non- vanishing terms of (4.6), i.e. only the terms of nt + rir + ne + — n < 0. Necessary 
formulae are developed in Ref. [10], and we have the mode decomposition of the singular self- force field schematically 
with vanishing ambiguous terms as 

. (4.8) 

I 

(Because of vanishing ambiguous terms, LHS and RHS are not exactly equal.) Again, the expansion coefficients are 
not continuous at the coincidence limit (f,r) {to,ro)- 

We subtract the singular self-force field from the bare self-force field. Because the singular behavior of these fields 
cancel each other, the sum of the harmonics becomes finite even after taking the coincidence limit, thus, the self-force 
is derived by 

I 

This shows that all we need in the end are the expansion coefficients at the coincidence limit {t,r) {to,ro ± 0), 
where the limit must be taken consistently. 



B. Radiation Reaction Formula 



The primary quc^stion to motivate the radiation reaction formula is a relation between the self-force and the balance 
formula. As we argue the self-force calculation in Subsec.IV A, a complicated regularization calculation is necessary to 
derive the self-force, while the balance formula is simply formulated without regularization. We also have a question 
which part of the self-force is described in the balance formula. We find the relation of the self-force and the balance 
formula in Ref. [13], and review some basic argument here. We also refer Ref. [14] for a hint to answer these questions, 
where the energy and angular momentum loss of the balance formula are derived from the radiative part of the 
self-force by taking an infinite time average. 

For preparation to discuss the radiation reaction formula, we first argue some symmetry property of a Kerr geometry 
and geodesies. Using constants of motion £"■, a = {E, Lz,C}, the geodesic equation in the Boyer-Lindquist coordinates 
becomes 

drV - /q!6'^^ 



-x) =^(^"'^)' [Tx) =^^'^-^'^^ ^'-''^ 

^ = Tr{£'^;r)+Te{S'^;0), g = $,(f ^ r) +$,(£«; 0) . (4.11) 

As we are interested in the inspiral stage of the binary, we only consider the case that r-motion is bounded, then, the 
r- and ^-motions become periodic [15] and we can expand them by discrete Fourier series as 

= ^ ^Hn) exp[znx6] , Xb = 27rO(,(A + A^) , (4.12) 

n 

where b = {r, 0} and r = z\ = z*^ , and, and fife are functions of £"'. Here we have two integral constants 

A^, but, since we can freely choose the zero point of A, only A'" — A^ can specify the geodesic. Because of these 
periodicities, one can freely add 1/51;, to A^. Suppose that the ratio of fJj, to fig is irrational, one can set A"" — A^ as 
small as possible, then, by an appropriate choice of the zero point of A, one can have A'' = A^ = 0. t- and </)-motions 
are simply integrated as 

z^ = Kc + Y, ^fc^"^ exp[inxfe] , Kc = < > A + C% (4.13) 

b,n 
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where c = {t, (j)} and t ^ z*, (f> = z"^, and, < Z'^ > and Z^^"-* arc function of £°. Here we have two integral constants 



C which specify the geodesic. Below, we call by primary constants, and, A and by secondary constants of 
geodesies. The Kerr geometry has the t and (^translation symmetries and, by these properties, one can set = 0. 
We find an interesting symmetry of the Kerr geometry. The metric is invariant under 

t^-t,r^r,e^e,(i)^-(i). (4.14) 

By this symmetry, the orbital parameter and the primary and secondary constants of a geodesic transform as 

A^-A, f-'^f", A''^-A^C"=^-C^ (4.15) 

Since the secondary constants can be set zero, the geodesic is invariant under this symmetry transformation and we 
call this by the geodesic preserving symmetry. 

We next discuss the self-force. Because the self-force is gauge dependent, one has to specify a gauge condition for 
this discussion. Here we restrict a class of gauge conditions such that the metric perturbation is derived as 

V(^)= j dx'G^,^.,.{t-t',cl>-(l>'-r,r',6,e')T^''-\x'), (4.16) 

where the Green function has an appropriate fall-off condition at |i — |r — r'| ^ oo. We call this by the physically 
reasonable class of gauge conditions because the Green function is invariant by the t and ^-translation symmetry of 
the Kerr geometry, and because one can naturally read out a spectral information of asymptotic gravitational waves 
in this class of gauge conditions. Instead of the coordinate components of the self-force, we consider the evolution 
of the "primary constants" by the self- force. The primary constants of a geodesic are defined by the Killing vectors 

r]a^^ and tensor ?7^^_j as i? = i]^v°',L = rj^v"^, C = r]^i^v"v^ /2, and we consider the self-force acting on the "primary 

constants" F°- = d£"'/dX derived as 

We note that (4.10) and (4.11) rely only on the definitions of whether they are constants or not, thus, (4.10), 
(4.11) and (4.17) describe the orbital evolution by the self-force. Because the self-force is induced by a geodesic to the 
leading order, it is a function of the primary and secondary constants and the orbital parameters. We also find that the 
self-force does not depend on C because of the t and (^translation symmetries, thus, we have F° = F°'{£°-,X^, C^; A). 
Under the gauge condition defined by (4.16), we find the self-force could be expanded by discrete Fourier series as 

F« = ^ £-(^r,ne) expltrirXr + mexs] • (4.18) 

nr,ne 

We consider to apply the geodesic preserving symmetry transformation (4.14) to the self- force. Because the trans- 
formation changes the direction of time, it also changes the boundary condition to derive the self-force. We derive 
the key identity, 

p{ret)a^^a^^b. ^) ^ _ p{adv)a ^^a ^ _ . ^ ^^_^g^ 

where we note the retarded boundary condition and the advanced boundary condition by {ret) and (adv) respec- 
tively to indicate the boundary condition for the self-force. Using the general form of the self- force (4.18), we have 
£(ret)a{n,,ng) ^ _g(adv)a(-n, -no) ^ Supposc we Calculate the self-force by the radiative (half-retarded-minus-half- 
advanced) Green function, we would have 

}^ ^g{ret)a{nr,ns) _^ g{ret)a{—nr,—ne)^ _ £{rad)a{nr,ne) (4 20) 

where we note the radiative boundary condition by {rad) to indicate the boundary condition for the self-force. This 
shows that we obtain the half of the self-force by using the radiative metric perturbation. The balance formula by 
Ref. [3] and Rcf. [14] discusses the infinite time averaged loss of the energy and angular momentum, corresponding to 
^(ret)a(o,o) _ ^(rod)a(o,o) ^ ^^^^ ^j^g prcscut derivation also includes the Carter constant. 

This result suggests why the balance formula was derived without a complicated regularization calculation as was 
required in the self-force calculation. Because the radiative Green function is defined as G''"'^ = (C^* — G°''^'")/2 
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[14], the singular parts of the retarded Green function and the advanced Green function cancel each other, thus, the 
regularization calculation is naturally done in the balance formula [3] and the radiation reaction formula [13]. 

Using the general form of the self-force (4.18), one can derive a general property of the orbit with the radiation 
reaction effect by integrating (4.10), (4.11) and (4.17). Here, it is important to observe that the self-force is derived 
by the linear metric perturbation. The divergence of the linearized Einstein tensor is known to vanish algebraically, 
thus, the stress-energy tensor of its source term conserves in the background. For this reason, the use of a linear 
metric perturbation as an approximation of the geometry is valid only when the deviation from a geodesic is small 
enough, and, in Ref. [13], we consider the orbit by a perturbation from a geodesic. By this perturbative analysis of 
the orbit, we find the evolution of the "primary constants" and the "secondary constants" by the gauge condition 
(4.16) can be written as 

= SS+<£''>X+ J2 f "^""'"^^ exp[m,xr + ineXe] , (4.21) 

rir.ng 

X>> = A(^+ < > AV2 + (A''("'-'"''U + A''^"-"^)) exp[zn^xr + inexe] , (4.22) 

nr,ne 

= <C^> AV2 + Yl (C'^^"'-"«U + C^("-"«)) exp[m,Xr + ineXe] , (4.23) 

where £q,\q, Cq are the initial values of the "constants". It is notable that < f >, < A*" > and < (7^ > determine 
the dominant contribution of (4.21), (4.22) and (4.23), which are derived only by the radiative metric perturbation. 
We also prove that these quantities are gauge invariant in the physically reasonable class of gauge conditions. 

Because the radiative Green function is a homogeneous Green function, there is a convenient method to calculate it 
even in a Kerr black hole [9] in the so-called radiation gauge conditions. A necessary calculation here is essentially a 
minor extension of that used for the balance formula, and we consider that a numerical method is already established 
[2]. 

C. Self-Force Calculation vs Radiation Reaction Formula 

Here we try to compare these two approaches to calculate the orbit. There could be two issues for comparison; 

1) Theoretical issue; whether the resulting orbit is physically acceptable or not 

2) Practical issue; whether the approach is practically available for generating templates for LISA project 

At this moment, we may say that both the self- force calculation and the radiation reaction formula still have the 
theoretical issue. 

Gravitational radiation reaction is a physically real object. Because the momentum flux carried by gravitational 
waves is well defined at the asymptotic flat region of the background spacetime, it is reasonable to expect the orbital 
evolution must be consistent with this effect. There is a wide-spread belief that the self- force includes the effect 
of gravitational radiation reaction. Although this looks a reasonable conclusion, we find that the relation between 
the gravitational self-force and gravitational radiation reaction in a usual metric perturbation scheme is not trivial 
because of the gauge freedom. As we shall argue in Sec.V, the self-force totally depends on a gauge condition in 
the whole time interval where a usual metric perturbation scheme is valid, thus, the gravitational self-force does not 
necessarily describe the effect of gravitational radiation reaction to the orbit as we have expected. In Sec. VI, we shall 
discuss that this contradicting problem appears because the usual metric perturbation scheme is not appropriate to 
consider radiation reaction to the orbital evolution, and we shall propose an adiabatic approximation of the metric 
perturbation where the physically reasonable class of gauge conditions is adopted. Unless the self-force calculation 
is done in the physically reasonable class of gauge conditions, we consider it necessary to argue that the self- force 
calculation in a gauge condition describe an orbital evolution by radiation reaction in an appropriate manner. 

On the other hand, the radiation reaction formula is derived consistently with the balance formula by using the 
physically reasonable class of gauge conditions, however, the radiation reaction formula gives us only the infinite 
time averaged part of the self- force acting on the "primary constants". Here we have a question. A general form 
of the self-force (4.18) shows that the self-force has an oscillating component other than the infinite time averaged 
component. It is not clear how this oscillating component contributes to the orbital evolution. (4.21), (4.22) and 
(4.23) show that the infinite time averaged component dominantly determines the orbital evolution, however, this 
analysis is based on a usual metric perturbation scheme. As we shall argue in Sec.V, the usual metric perturbation 
is not appropriate to consider the orbital evolution, and we shall propose an adiabatic approximation of the metric 
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perturbation in Sec. VI. Then the question remains on a possibility that the oscillating component might be important 
to consider the orbital evolution. 

As we shall discuss in Sec. VI, we consider the radiation reaction formalism under the adiabatic approximation of 
the metric perturbation. We find that the infinite time averaged part of the self-force actually determines the orbit 
during the time scale where the adiabatic approximation of the metric perturbation holds. We also find that the 
oscillating part of the self- force can be eliminated by a special choice of a gauge condition in Ref. [16]. 

As for the practical issue, still the self-force calculation has a lot to investigate. There are a number of studies of a 
linear metric perturbation induced by a point particle [2]. however, most of these studies calculate gravitational waves 
at infinity where both a numerical and a semi-analytic methods arc established. On the other hand, the self-force 
calculation needs a mode decomposition of a linear metric perturbation along the orbit in order to derive the bare 
self-force field (4.5) and a reliable technique for this calculation has not been demonstrated in general. 

If the background is a Schwarzschild black hole, one can use the Zerilli-Regge- Wheeler formalism [17] to derive a 
mode decomposition of the bare self-force field for a general orbit in principle, and we have a result for a circular 
orbit or a radial orbit, not for a general orbit. But, when the background is a Kerr black hole, a method to calculate 
a metric perturbation is just proposed [8], and the self-force calculation with this idea is not yet formulated. A more 
crucial problem here is that it is not clear whether the self-force calculation using this metric perturbation would 
describe the effect of gravitational radiation reaction in a physically reasonable manner because the gauge condition 
along the orbit is not well understood. 

We also comment that the self-force calculation is a point splitting regularization, and that the singular self-force 
field might not be optimally subtracted from the bare self-force field. As a result, a numerical convergence of the 
self-force calculation might not be fast in general. 

On the other hand, the radiation reaction formula can be seen as a non-trivial extension of the balance formula. 
We only need a radiative metric perturbation and a convenient method for this calculation is known [9]. Various 
numerical codes of the balance formula arc already made successfully [2] for either a circular orbit or an equatorial 
orbit around a Kerr black hole, and we are coming to have a result for a general orbit, therefore, we consider this 
approach is practically promising. There is also known a semi-analytic technique which may substantially increase 
the efficiency of the code [18] . We also comment that we have the most optimal convergence in this approach. 

V. SELF-FORCE? A MYTH? 

So far, we frequently noted that the self-force does not necessarily include the effect of gravitational radiation 
reaction, which is obviously against our motivation to study the self-force problem. Because it has not been seriously 
considered whether the self-force includes the effect of gravitational radiation reaction, we shall devote this section to 
discuss this problem. 

Gravitational radiation reaction is physically real since we can define a momentum flux of gravitational waves at 

the asymptotic flat region of the background metric, however, we find that the self-force is entirely gauge dependent 
and that it could vanish by a special gauge choice along the orbit. It is well-known that the self-force could vanish 
at an orbital point by a gauge choice, but we considered this might not be a problem because the self-force would 
have non- vanishing components of gravitational radiation reaction by taking an infinite time average of the scJf-forcc. 
In Subsec.IVB, we discuss that this is the case in the physically reasonable class of gauge conditions and that non- 
vanishing components are proven to be gauge invariant in this gauge class. However, a usual metric perturbation 
scheme can describe a gravitational evolution of the system only for a finite time interval, and, in this whole time 
interval, the self- force could even vanish by a gauge choice along the orbit as we show in Subsec. V A. This suggests that 
the relation between the self- force and gravitational radiation reaction is entirely gauge dependent, and is not trivial 
as we expected before. Then we have questions; how the self-force is related with gravitational radiation reaction, 
and what gauge conditions we should take to calculate a self- force. We consider these questions in Subsec. V B 

A. Metric Perturbation 

We suppose that a regularization calculation for a self-force is formally possible by a matched asymptotic expansion 
as we discuss in Subsec. Ill C, and consider only a metric of the far-zone expansion. The metric of the far-zone 
expansion is a sum of a regular vacuum background metric and its metric perturbation induced by a point particle 
moving on the background metric. We suppose the point particle has an appropriate structure [7] so that the metric 
of the far-zone expansion consistently matches a metric of the near-zone expansion. 
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Suppose wc use a usual metric perturbation scheme for the calculation of the metric of the far-zone expansion. We 
expand the metric and the stress-energy tensor of the point particle with a small parameter m/ L as 



.9m. = g'iS^ + {^/L)h%l + {^/Lfh^!}} +■■■, (5.1) 
T^-^ = (m/L)r(i)^'' -h {m/LfT^'^'^i''' + , (5.2) 



where g^^v^ is the vacuum background metric. For a valid perturbation, we have 



0(1) > (m/L)/iW > {m/Lfh^^} >■■■, (5.3) 
0(1) > im/L)T^^^'"' > {m/LfT^^^'"' > ■■ ■ . (5.4) 



One can expand the Einstein equation with m/L, and we schematically have 



^ (5.5) 
G(i)a"'[/i(2)] + g(2)'^'^[/i(i), = t(2)'^'' , (5.6) 



where G^^^^''[/i] and G^2'^''[/i] are the terms linear and quadratic in h^i^ of the Einstein tensor C^^lg + h], respectively. 
As we discuss in Subsec.IIIC, T^'^^'^" is a usual stress-energy tensor of a monopole particle. It is important to 

note that the linearized Einstein equation is algebraically divergence free with respect to the background metric as 
QWfj-'^^ — r(i)^'^.^ = 0, thus, in this metric perturbation scheme, the particle as a source of h'l^,} must move along a 
geodesic in the background metric g^^i?' for a consistent solution of (5.5). An explicit form of 

may be derived 

by doing a matched asymptotic expansion of metrices of the far-zone expansion and the near-zone expansion to this 
order. But, at least, one may consider that it has a term of the monopole particle deviating from a geodesic because 
one can derive the MiSaTaQuWa self-force from LHS of (5.6) by a mass renormalization [1]. Thus, for the expansion 
of the stress-energy tensor, we must also consider the expansion of the orbit as 

z^'{T) = z^'x^^'^ir) + {m/L)z^^^i'{T) + • • • , (5.7) 

where r is defined to be a proper time in the background metric and is a geodesic of the background metric. 

For a valid perturbation, we have 

0(1) > (m/i)^(i)'' > • • • . (5.8) 

Intuitively speaking, the orbit deviates from a geodesic by gravitational radiation reaction, and eventually the condition 
(5.8) would be violated, then the metric perturbation in this perturbation scheme would fail to approximate the system. 

We consider a gauge transformation in the time scale where the metric perturbation still approximates the system. 
The gauge transformation is defined by a small coordinate transformation ^ = -\- (to/L)^'^(.e) and the orbit 
transforms as z^(t) z^{t) = z^(t') + (m/ L)S^'^{z{t')), where t' is related with the orbital parameter t so that t 
remains a proper time of the new orbit z^{t) in the background metric. Applying this to the perturbation expansion 
of the orbit (5.7), we find the gauge transformation of {m/L)z^^^'^ as 

(m/L)z(i)^(r) ^ (m/L)z(i)^(r) = (m/L)z(i)''(r) + v^^<^'^^'{t' - r) + (m/L)C^(z(''9) (r)) , (5.9) 

where u^'s)^ = dz^^^^^ /dr. We find that we could eliminate {m/ L)z^^^^{t) if the gauge transformation along the orbit 
satisfies 

(m/L)e^(z(''9)(T)) = -(m/L)z(i)^(r) +t;^5T, (5.10) 

with an arbitrary small function 6t ^ 0{m/L). There exists a gauge transformation which satisfies this condition for 
the whole time interval where the metric perturbation scheme is valid because it is smaller than 0(1), and the orbit 
becomes a geodesic of the background metric in this gauge condition. 

Because the self-force is an self-acceleration to deviate the orbit from a geodesic, {m/ L)z^^'>^{t) = in this gauge 
condition means that the self- force entirely vanishes for the whole time interval where the metric perturbation scheme 
is valid. This extreme example suggests that the sclf-forcc is a totally gauge dependent object in the usual metric 
perturbation scheme, and, even a time average of the self-force over a longest time scale where the metric perturbation 
is valid is gauge dependent in general. 
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The essential reason why the self-force could not include the effect of gravitational radiation reaction is that this 
perturbation scheme only allows a small deviation from a geodesic as in (5.8). Because of this, one could always 
eliminate the deviation by a gauge transformation as in (5.9). We note that this problem may not be solved by 
calculating a non-linear metric perturbation because the problem comes from the metric perturbation scheme. We 
consider a solution of this problem is to modify the metric perturbation scheme so that one can describe a non- 
perturbative orbital deviation from a geodesic. If the orbit can deviate from a geodesic non-perturbatively, one 
cannot eliminate it by a gauge transformation. 



B. Energetics of the Orbit and Radiation 

A hint to imderstand why such an unexpected thing happens can be seen in Rcf. [13]. In Sec. Ill of Rcf. [13], we 
argue that the orbital energy does not decrease monotonically by the emission of gravitational waves. We consider 
that gravitational radiation has its own energy and that the self-force describes the interaction between the orbital 
energy and the radiation energy. In fact, one can define an effective stress-energy tensor of the orbit and gravitational 
radiation [1]. As in the previous subsection, we only consider the metric of the far-zone expansion, and we suppose 
that the metric is a sum of the vacuum background metric and its perturbation induced by a point particle with an 
appropriate internal structure [7] as 9^,^ = .g^'i?' + hfj_^. 

The Einstein equation can formally written as G^i^[/i] + G^l+'[/i] = T^,,, where C'j!^ and o'j^/^ are the linear 
terms and the rest of the Einstein tensor with respect to /i^^ , and T^i, is the stress-energy tensor of the point particle. 

Because is algebraically divergence free, one can define a conserved stress-energy tensor in the background metric 
by 

T^"' = T'''' ~G^^+^^"'[h]. (5.11) 

We consider that the first term and the second term of (5.11) represent the stress-energy tensors of the particle and 
gravitational radiation respectively, and each of these are not conserved in the background by themselves. We suppose 
that the background metric g^t^'' is a Kerr black hole and has a timelike Killing vector then one can define a total, 
orbital and radiation energy in the background metric as 

Eitot) ^ £;(or6) ^ ^(rad) ^ ^(orb) ^ _ j rf^^^^T^/^^ ^ ^ j dY,i,^^G^^+^^'' , (5.12) 

where the surface integration is taken over a spacelike hypersurface bounded by the future horizon and the future null 
infinity of the background black hole metric. We note that the radiation energy is not necessarily positive. 

By integrating (5.11) over a small world tube surface around the orbit, one can derive the MiSaTaQuWa self-force 
[1], which shows that the self- force describes the interaction between the orbital energy and the gravitational radiation 
energy. A gauge condition for the self-force could be interpreted as a small arbitrariness in defining the orbital energy 
and the radiation energy in (5.12) [16]. In order to see this, we apply the metric perturbation scheme of Subsec.VA 
to (5.12), and we have 

= (m/L)£(°'''')(i) + {m/LfE^'"-^^^^^ + ■■■ , 
{m/L)E^"'-^^^^^ = - J (n:^^^{m/L)T^^^i"' , {m/LfE^"^'''^^^'^ = - J d'E^^,{m/LfT^'^^'"' , • • • (5.13) 

Eirad) ^ (to/L)2^(°'-'')(2) + . . . , {m/ if E^'-"^^^'^^ = J dE^^^m/ L)^ G'^^^^"' [h^^\ h^^^] , • • • , (5.14) 

In the usual perturbation scheme, (m/L)T(-'^)'"^ conserves by itself, as a result, becomes a constant of motion. 

By a gauge transformation around the orbit, one can show that the orbital energy and the radiation energy transform 
as 

(m/L)2^(°'''')(2) ^ (m/L)2^(°'-'')(2) = {m/LfE^'"'''^^^^ + {m/LfSE, (5.15) 

(m/L)2£(™'*)(2) ^ (m/L)2£(™<^)(2) = {m/LfE^'''"^'^^'^^ - {m/LfSE , (5.16) 

where 5E is arbitrary, and the orbital energy is entirely gauge dependent by itself as we discuss in the previous section. 
It is important to note that is gauge invariant-'- [16]. 



Because of this property, one can use this as a criteria for a last stable orbit of an extreme mass ratio binary. 
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Although wc do not yet propose a new metric perturbation scheme, we may extend the argument of the orbital energy 
and the radiation energy since they arc defined in a non-perturbative manner as (5.12). The balance formula is derived 
by gravitational waves at the future null infinity and the future horizon of the background black hole geometry, and it 
describes radiation reaction to the total energy rather than the orbital energy [16]. In the usual metric perturbation 
scheme, one can decrease only the radiation energy by radiation reaction while keeping the orbital energy constant by 
the gauge freedom. Although this is mathematically allowed in the usual metric perturbation scheme, if the system 
continues to lose the total energy by radiation reaction, it will eventually be difficult to understand the orbit in a 
physically reasonable way. Since we keep the orbital energy constant, the radiation energy will inevitably decrease 
substantially by radiation reaction, as a result, the amplitude of the metric perturbation would be non-perturbatively 
large and one cannot use the coordinate system of the background metric as a reference to observe the orbit. 

In this case, only numerical relativity may be able to describe the evolution of the system, and, even if it is possible, 
it may be difficult to understand the system. Thus, for this technical advantage, we consider it reasonable to keep the 
amplitude of the metric perturbation small and we should decrease the orbital energy by radiation reaction rather 
than the radiation energy. We consider that this would be a reasonable least criteria for a convenient gauge condition 
to calculate the self-force. Although we do not have a simpler criteria for a convenient gauge condition, the physically 
reasonable class of gauge conditions (4.16) satisfies this criteria. Because the self-force in this gauge class is consistent 
with the balance formula, only the orbital energy decreases by radiation reaction and the radiation energy is just 
oscillating, as a result, the metric perturbation would not grow non-perturbatively large. 

Although the argument of energetic shows us an importance of the self-force, this is just a qualitative argument. 
In the next section, we propose a new perturbation scheme and we formulate a method to calculate a gravitational 
evolution of an extreme mass ratio binary in this picture. 



In this section, wc review our recent progress on the radiation reaction formalism in Ref. [16]. Wc discuss in 
Sec.V that the self-force does not necessarily include the effect of gravitational radiation reaction because of the 
perturbation scheme. In order to break this limitation, we introduce an adiabatic approximation to this problem. 
The adiabatic approximation is well-known in classical mechanics, however, the application to gauge field theory such 
as a gravitational perturbation is not so common. 



From the result of the radiation reaction formula, the orbital deviation becomes 0(1) when A ~ 0{{m/L)~^^^). 
This time scale is called a dephasing time of a orbit. The usual metric perturbation scheme is valid only within this 
time scale as wc discuss in Sec.V. Our purpose here is to extend this time scale of validity by modifying the metric 
perturbation scheme. For an explicit discussion, we use the physically reasonable class of gauge condition (4.16). 

We denote the orbital constants of a geodesic by 7 = {£"■, X^, C}. Because a linear metric perturbation is induced 
by a geodesic, it is a fimction of 7 as hf^i,{x) = h^uAj', x). In the radiation reaction formula [13], we use an adiabatic 
approximation to the orbit. In this approximation, we approximate the orbit by a geodesic at each instant, and 
consider the evolution of the "orbital constants" by the effect of gravitational radiation reaction, thus, we may write 
the "orbital constants" as functions of the orbital parameter as 7(A). We consider to extend this idea to the metric 
perturbation. We foliate the spacetime by spacelike hypersurfaces which intersect with the orbit. We define the 
foliation function f{x) by the orbital parameter at the intersection of the orbit and the surfaces as f{z{\)) = A. We 
define an adiabatic linear metric perturbation by using the linear metric perturbation on the foliation surface induced 
by the geodesic 7(/) as 



Here we do not specify an explicit form of the foliation function, as a result, an adiabatic linear metric perturbation is 
not defined uniquely. We consider that we may need a constraint of the foliation function for an adiabatic non-linear 
metric perturbation, however, it docs not change the result of the following discussion. We consider that the adiabatic 
linear metric perturbation is a non-trivial extension of the linear metric perturbation. 

Under the physically reasonable class of gauge conditions, a formal expression of the tensor Green function becomes 



VI. ADIABATIC EXTENSION 



A. Adiabatic Metric Perturbation 



(6.1) 



G^.^>Ax^x') = 9t."Mr, 9; r' , 9') exp[-za;(i - t') + *m(</. - 4>')] . 



(6.2) 
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This gives a general form of the hnear metric perturbation in this class of gauge conditions induced by a geodesic 7 

as 

V(7;^)= E fc^"''"'"'-'"^H7)/i|r.'"'"'-'"^H^^^,^)exp[-ia;t + *m</.], (6.3) 

(7) = J dAfc('"'™'"'-'"^) (f exp[iuJKt - imn^ - iUrXr - inexe] ■ (6.4) 

We construct the adiabatic linear metric perturbation (6.1) with the general form of the linear metric perturbation 
(6.3). When we operate the linearized Einstein operator, wc have 

G« n = V[7(/)] + AW [h'^'^] , (6.5) 

where an extra term A^"^j appears because the adiabatic linear metric perturbation is not induced by a geodesic of 
the background. It is notable that T'^'^['y{f)] is a stress-energy tensor of a point particle moving along the orbit with 
the effect of gravitational radiation reaction. 

The adiabatic metric linear perturbation solves the Einstein equation to an accuracy of 0(m/L) as long as 0{m/ L) > 

A^i^ holds. Since A^^} is proportional to the A-derivatives of the orbital "constants" 7(A), the validity of the adiabatic 
metric perturbation depends on how the orbit evolves by a self-force. In the next section, we discuss the orbital 
evolution under the physically reasonable class of gauge conditions. Wc find that the adiabatic metric perturbation is 
a well approximated solution of the linearized Einstein equation in the radiation reaction time scale 0((m/L)~^) > A, 
hence, the time scale of validity is much longer than that of the usual metric perturbation scheme. 



B. Adiabatic Evolution of the Orbit 



Different from the usual perturbation scheme one can not use the perturbation of the orbit any more, and we need 
to calculate the orbit in a non-perturbative manner. Because of a technicality of a non-perturbative calculation, we 
only summarize the result and we suggest the readers to refer Ref. [16] for the details. 

The self-force defined by the adiabatic linear metric perturbation is simply derived as 

d_^a ^ pa^^a^^^^ ^'^^A); A) + 0{{m/Lf) , (6.6) 

and we can use the result of the original radiation reaction formula [13] in the radiation reaction time scale A < 
0{{m/L)~^). As for the "primary constants", one can still deal with the evolution by a perturbation, and we find 
that 5"^"'°) dominantly determines the orbital evolution. The subdominant parts of the "primary constants" are 
0{m/L) and we may not see these effect by gravitational waves. The evolution of the "secondary constants" becomes 
non-perturbative beyond the dephasing time, however, we find the perturbation results of (4.22) and (4.23) are 
qualitatively correct. The dominant part of their evolution is described only by £'«(o>o)^ ^nd the subdominant part 
grows linearly in A. 

In summary, we find the qualitative behavior of the orbital evolution as 

^£-=0{m/L), ^X' = 0{{m/m, j^C^ = 0{{m/ L)X) . (6.7) 

This gives us the qualitative estimate of (6.5) and we find A^i^ = 0{{m/L)'^X). Hence, the adiabatic linear metric 
perturbation is valid in the radiation reaction time scale 0{{m/L)~^) > X. 



C. Radiation Reaction Gauge 

The validity of the adiabatic linear metric perturbation depends on a behavior of A\}} in (6.5). Because A\}} is gauge 
dependent, we may consider to extend the validity by using a remaining gauge freedom in the physically reasonable 
class of gauge conditions. 

Using the adiabatic linear metric perturbation, we find that the "primary constants" still evolve perturbatively as 

= £S+ > X+ f "^"^'"'^ exp[m^xr + ineXe] ■ (6.8) 
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By a gauge transformation — > a;** + the "primary constants" transforms as f " ^ f " + SS'^ where 

SS^/^ = -vTv^Q + r^S^/^r , (6.9) 
^^"^ = -V^yV'^^^ . (6.10) 
In Ref. [16], we find that, by an appropriate gauge choice, one can eUminate the oscillating part of (6.8) as 

= + s£S)+ < f > A + oi^ih) . (6.11) 

Here, one has to introduce a small shift of the initial values 6£q so that the gauge transformation behaves as ^ = 
0{m/L). In this gauge condition, the self-force has only the radiative part as 

F»(£:"(A), A^(A), C^(A); A) = < f " > (£:"(A)) . (6.12) 

For this reason, we call this by a radiation reaction gauge. 

The self-force is expected to have a conservative part as well as the radiative part, however, in this gauge condition, 
the conservative part is integrated out to be a rcnormalization of the initial values £q ^ f g + 6£q. We also note 
that a radiation reaction gauge condition is applicable to an orbit of a spinning particle. It is known that the orbit of 
the spinning test particle deviates from a geodesic [19] by a coupling of its spin and the curvature of the background 
metric. This effect can also be renormalized to the initial values in the radiation reaction gauge condition. 

Using a post-Newtonian estimation [20], the self-force can be estimated to be 0{m/L x v^), where v'^ ~ 0.1 around 
a last stable orbit. This is smaller than that in a general gauge condition of the physically reasonable class. As a 
result, one can predict the orbit in the time scale A < 0{{'m/Lv^)~^), which may correspond to several years for a 
promising target of LISA. 



VII. SUMMARY AND FUTURE PROSPECT 



It is widely believed that a self-force could describe an orbital evolution with the effect of gravitational radiation 
reaction. For this reason, the so-called Capra community were trying to develop a method to explicitly calculate 
the self-force acting on a particle orbiting around a Kerr black hole, and we arc having a success in developing a 
regularization calculation method either by the self- force calculation [10] or by the radiation reaction formula [13]. 
On the other hand, we come to have a concern on what physical information the self-force actually carries because 
the self-force is gauge dependent. For this concern, it was discussed that it might be necessary to develop a non-linear 
metric perturbation to grasp a physical meaning of the self-force. 

A second order metric perturbation may be derived consistently only when we consider an orbit with the effect 
of gravitational radiation reaction, and it gives us a gauge invariant information by an asymptotic gravitational 
waveform, for example, by a modulation of a gravitational wave phase. We consider that we may see the radiation 
reaction effect from the waveform, however, it has nothing to do with the radiation reaction effect to the orbit because 
radiation reaction to the orbit is purely gauge dependent. We find that one can define a radiation energy, an orbital 
energy, and a total energy as a sum of the radiation and orbital energy. Radiation reaction we usually consider (say, 
by the balance formula) is a reaction acting on the total energy, but, the orbital motion is related to the orbital energy. 

In Sec.V, we discuss that the gauge freedom is interpreted as a small ambiguity in separating the total energy into the 
radiation energy and the orbital energy. This separation is entirely arbitrary in the usual metric perturbation scheme, 
however, we discuss that, without some constraint on this separation, we may have a difficulty in interpreting the 
orbit. When we continue the evolution beyond the usual metric perturbation scheme, the metric perturbation would 
grow non-perturbatively large and, as a result, the background metric cannot be used as an approximated reference 
to track the orbit. We argue that the physically reasonable class of gauge conditions is, at least, a reasonable gauge 
choice to avoid such a non-perturbative situation in Sec.V, and we show that an adiabatic extension of the linear 
metric perturbation is possible in Sec. VI. 

A possible question here is whether there is a wider class of gauge conditions which allows an adiabatic extension 
of a linear metric perturbation. Especially, we introduce the physically reasonable class of gauge conditions because a 
linear metric perturbation is derived by a Fourier-Harmonic decomposition, and it is not clear whether an adiabatic 
extension of a linear metric perturbation by a time domain calculation is possible or not. If it is not possible, the 
self-force calculation by the time domain calculation may not predict an orbital evolution in the radiation reaction 
time scale. 

Different from the self-force calculation, the radiation reaction formula conveniently calculates only an infinite time 
averaged part of the self-force. We find that this part is actually enough for the prediction of an orbit in the radiation 
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reaction time scale. We find the rest of the self-force only makes a small change of the orbit, which can be totally 
eliminated by a choice of the gauge. 

A crucial problem to calculate gravitational waveforms for the LISA project is whether the radiation reaction time 
scale is sufficiently longer than an observation time of an astrophysically expected target. If we calculate the self-force 
in the physically reasonable class of gauge conditions, the radiation reaction time scale would around several months 
in general. Since an observation time of LISA would be around several years, there would be a case that one cannot 
calculate reliable waveforms. We find that one could modify the radiation reaction time by using the remaining gauge 
freedom in this class. We propose the so-called radiation reaction gauge where the self-force has only the radiation 
reaction component, which is derived by the radiation reaction formula. In this gauge condition, we find that we 
would have reliable waveforms of several years, which is sufficient for the present LISA project. 

Although we consider the radiation reaction formula with an adiabatic linear metric perturbation may be sufficient 
for LISA project, it remains a great theoretical challenge to calculate an orbit longer than the radiation reaction time 
scale. 
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